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Abstract 


We derive the quantum analogue of a Backlund transformation for 
the quantised Ablowitz-Ladik chain, a space discretisation of the nonlin¬ 
ear Schrodinger equation. The quantisation of the Ablowitz-Ladik chain 
leads to the g-boson model. Using a previous construction of Baxter’s 
Q-operator for this model by the author, a set of functional relations is 
obtained which matches the relations of a one-variable classical Backlund 
transform to all orders in h. We construct also a second Q-operator and 
show that it is closely related to the inverse of the first. The multi- 
Backlund transforms generated from the Q-operator define the fusion ma¬ 
trices of a 2D TQFT and we derive a linear system for the solution to the 
quantum Backlund relations in terms of the TQFT fusion coefficients. 

1 Introduction 

In the context of classical integrable systems the main interest in the construc¬ 
tion of Backlund and Darboux transformations is their application in the con¬ 
struction of soliton solutions [13 [20]. Compared to the classical Backlund trans¬ 
formations the discussion of their quantum cousins started more recently. In 
1992 Gaudin and Pasquier [H] constructed for the quantum Toda chain an ana¬ 
logue of Baxter’s Q-operator [3 Ch. 9-10] and showed that in the semi-classical 
limit /i —>■ 0 the similarity transformation O —^ Q{u)OQ{u)~^ is a Backlund 
transform of the Toda chain. An introductory account to this and further re¬ 
sults can be found in Sklyanin’s lecture notes p7] . 

1.1 The Ablowitz-Ladik chain and its quantisation 

In this article we are interested in the quantum analogue of Backlund transforms 
for another integrable system: the Ablowitz-Ladik chain [T], 
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Figure 1: Vertex configurations for the T-operator. Consider a square lat¬ 
tice where we assign a nonnegative integer to each lattice edge, the statistical 
variable. At each vertex we only allow for particular configurations, for the T- 
operator the horizontal edges can only take values 0 or 1, below are the vertex 
configurations in terms of non-intersecting paths. 


where we consider periodic boundary conditions '4’j+n = V’j and '4’j+n = 

The equations are a space discretisation of the following system of coupled 
PDFs 

f dtip = d'^ip — 

\ dt^* = -dlip* + ■ 

After changing to imaginary time, t ^ it with i = \/—l, this system of PDFs 
allows for a reduction, ip* = ±'0 with ip denoting the complex conjugate of ip, to 
the nonlinear Schrddinger (NFS) equation, —idtip = d'^ip^2ip\ip\^. Backlund and 
Darboux transformations for the system (HID and its reduction to the discrete 
NFS system can be found in e.g. [HllilllHllSl] and references in loc. cit. 

Kulish considered in m a particular quantisation of the Ablowitz-Fadik 
chain leading to the g-boson model with Hamiltonian 

n 

H = - - 2(1 - ( 1 - 2 ) 

i=i 

see also mw as well as references therein. Here the /3j , /3* ’s are the generators 
of a g-deformed version of the oscillator or Heisenberg algebra with q being the 
quantisation parameter (see the definition (I3.2|l in the text) and Nj is the particle 
number operator at site j. 

In the Fock space representation and with periodic boundary conditions on 
the lattice, Pj+n = the model (jl.2ll can be solved via the 

quantum inverse scattering method [9]. Analogous to the classical case, one 
defines a quantum Fax operator L{u) and the Hamiltonian (II.2|) can then be 
understood as a particular element in a commutative algebra generated from the 
commuting transfer matrices T(u) = Tr L(u) of an exactly solvable lattice model 
in the sense of Baxter [2]; see Figure [T] for the vertex configurations defining the 
model. 
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Figure 2: Vertex configurations for the Q^-operators. Consider a square lat¬ 
tice where we assign a nonnegative integer to each lattice edge, the statistical 
variable. At each vertex we only allow for particular configurations, for the 
Q^-operator (bottom left) we require that b > c and a + b = c + d. For Q~ 
(bottom right) we impose instead the conditions a -I- 6 > c and a + b = c + d. 
On the right of each vertex is an interpretation of these conditions in terms of 
non-intersecting lattice paths. 


1.2 Baxter’s Q-operators and quantum Backlund maps 

In this article we discuss two analogues of Baxter’s Q-operator for the 
quantised Ablowitz-Ladik chain (O with periodic boundary conditions. The 
construction of these Q^-operators corresponds to the dehnition of two addi¬ 
tional exactly solvable statistical mechanics models defined on a square lattice; 
the corresponding vertex configurations are depicted in Figure [2] We show that 
the transfer matrices T{u), Q^(u) of all three models are related via functional 
relations, Baxter’s famous TQ-equation and a quantum Wronskian relation for 
Q^. The Q^-operator has been constructed earlier by the author, the construc¬ 
tion of Q~ and the resulting functional relations (|4.8I1 . (14.1311 in the text are 
new. We show that in the limit n —>■ oo the second solution Q~ becomes the 
inverse of . 

Employing the Q^-operator we consider the associated similarity transfor¬ 
mations of the g-boson fields, 

l3^^]3j{v)=Q+{v)l3fi+{v)~^ and /3* ^*(u) = Q+(z;)/3*Q+(?;)"^ . 

The first main result of this article is the proof that the transformed fields (3^ (v) 
and /3j(v) obey the functional relations (c.f. Theorem 15.21 in the text) 






and 



- 1 ), 


(1.3) 


which allow one to compute ]3j{v) and Pj{v) recursively via a power series ex¬ 
pansion in V. This result is physically significant, since this Backlund map (and 
the one induced by the adjoint of the (^'''-operator) describes for small time 
steps 0 < n = At <C 1 the discrete time evolution of the quantum Ablowitz- 
Ladik chain. Surprisingly, the quantum relations (El match exactly (i.e. to all 
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orders in h) the relations of the classical fields '0^, under the analogous clas¬ 
sical one-variable Backlund map B^{v) : >->■ '0^ (z;)) considered 

by Suris in [28] for the system (HD. 

In comparison, the match between Baxter’s Q-operator and the classical 

Backlund map for the Toda chain is established by formulating Q as an integral 

operator and then identifying its integral kernel in the limit h ^ 0 with the 

generating function of the Backlund transform; see [24] , Our approach avoids 

the generating function altogether and directly arrives at the functional relation 

m which allows one to compute the image of the quantum transform, the 
^ ____ 
quantum variables jSj, Pj, via recurrence; see Section [oTT] 

1.3 Multivariate Backlund maps and 2D TQFT 

Once a one-variable Backlund map is constructed one can consider multivariate 
transforms via composition, {xi)oB^{x 2 )o- ■ ■oB^{xn-i)- These multivariate 
transforms are central to Sklyanin’s separation of variables approach; see e.g. 
[18] and, for a discussion of the quantum separation of variables approach, [22] 
as well as references therein. Here we wish to highlight a novel aspect of these 
multivariate transforms: they generate the fusion matrices of a 2D topological 
quantum field theory (TQFT) constructed in [15) . 

Consider matrix elements of a product of Q'^'-operators which is the quantum 
analogue of the above multivariate Backlund transform, 

n—1 

(A| n Q^ixi)\n) = ... ,Xn-i-,q) , (1.4) 

where A, /r, v are partitions labelling particle configurations of the g-boson 
model (see Eqn (13.1511 in the text) and the P\s are a special basis in the ring 
of symmetric functions, the so-called g-Whittaker functions which are a special 
case of Macdonald’s functions m Ch. VI]. Since the Q+-operators for different 
xPs commute, reflecting the analogous property of the classical Backlund trans¬ 
forms, the above expansion is well-defined. The expansion coefficients iV^j^(g) 
are polynomials in q with integer coefficients and define the fusion in a 2D 
TQFT, i.e. they are the values of the pair of pants cobordism shown in Figure 
|3| which can be seen as a 2D analogue of a Feynman diagram describing the 
fusion of particles in a QFT. The TQFT defined via (II.4p is a q-deformation of 
the S't/(n)-WZNW fusion ring, that is the constant term in is 

given by the operator product expansion of primary fields in WZNW conformal 
field theory. 

The second main result of this article is that the matrix elements of the 

~ ~ jjt 

Backlund transformed g-boson fields Pj{v) and Pj{v) can be computed in terms 
of the TQFT fusion coefficients, thus relating the quantum Backlund transform 
of the Ablowitz-Ladik chain to fusion in a 2D TQFT; see Prop 16.41 in the text. 
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Figure 3: Depiction of the ‘pair-of-pants’ cobordism which fixes the fusion prod¬ 
uct in a 2D TQFT. The dotted lines inside are the corresponding conventional 
Feynman diagram. 

1.4 Outline of the article 

The outline of the article is as follows. In Section 2 we review some aspects of 
the classical Ablowitz-Ladik chain, in particular we recall a one-variable family 
of Backlund transforms considered by Suris in |28] which is quantised in terms 
of the Q^-operator in Section 4. Before that Section 3 reviews the necessary 
background on the quantisation of the Ablowitz-Ladik chain, the g-boson al¬ 
gebra and the algebraic Bethe ansatz solution of the g-boson model. Section 
4 gives the construction of the Q^-operators and the derivation of Baxter’s 
TQ-equation as well as the resulting quantum Wronskian relation. Section 5 
contains the main result, the discussion of the quantum Backlund transform 
and the derivation of the functional relations (IE3. In Section 6 we then relate 
these relations to the 2D TQFT defined in terms of the Q-operator. Section 7 
contains a discussion of the results and outlook towards future work. 

2 The classical Ablowitz-Ladik chain 

Interpreting the discrete fields {ipj}j&^ and {'0j}jez„ as components of vectors 
in K^", the classical Ablowitz-Ladik chain allows for the following Poisson 
structure (see e.g. [17] and [28] 1 



Note that the discrete fields tpj, ijj* are not a pair of canonical variables. The 
time evolution can be described as Hamiltonian flows, = {H, ipj} and 
dt-tp* = where (assuming < 1) 



( 2 . 2 ) 


Initially, one can consider ipj and ip* as independent variables, but as explained 
in the introduction, the system reduces to the NLS model when changing to 
imaginary time, t —>■ it, and identifying ip* as the complex conjugate of ipj. 
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The integrals of motion of the system (HU are obtained from the spectral 
curve det{L{u) — X) — 0, where the Lax operator L{u) is given as a product of 
the following local Lax matrices, 


L{u)=Lr.{u)---L2{u)L,{u), L,{u)=(^^^ 

Here we have chosen not to take the standard form of the Lax matrices, see 
e.g. [HIlTlISl] and references therein, but instead match the conventions in m 
which turn out to be convenient to discuss the quantum case. There are two 
spectral invariants, the determinant 


detL{u) = (2.4) 

j=i 

and the trace 

n 

T{u) = Tr L{u) = u^Tr . (2.5) 

r—0 

Using the classical analogue of the Yang-Baxter equation one proves the follow¬ 
ing [12]: 

Proposition 2.1 {Tr, Tr'} = 0 for all r,r' = 0,1,... ,n. 

The integrals of motion Tr create higher flows and their existence implies 
Liouville integrability of the system (HID- The Hamiltonian p.2ll creating the 
physical time flow (HU is given by 

H =-Ti-Tn-i+2Tr{\nL{l)) (2.6) 

Thus, we note that the physical time flow separates into three different com¬ 
muting flows generated by Ti = Tn-i = j+i and the 

determinant (EU; see the discussion in [28]. For instance, the flow generated 
from —Ti gives rise to the differential-difference equations [28l Eqns (2.7) and 
( 2 . 8 )] 

dtip* = -fj^+iil - 

The classical Backlund transforms which we discuss in the next section can be 
understood as time-discretisation of this particular flow. 

2.1 A one-variable family of Backlund transforms 

In this article we are interested in one-variable families of Backlund transfor¬ 
mations B^{v) : {ip. Ip*) 1 -^ {ip, ^ ) which correspond to the flows discussed by 
Suris in [28l Prop 1 and Prop 4] with v playing the role of the discrete time pa¬ 
rameter vZ. We summarise here some of the results from loc. cit. for the case of 
periodic boundary conditions to keep this article self-contained and to introduce 
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our notation. The novel aspect in this article is that we relate these discrete 
time flows to the Q^-operator which we introduce in a subsequent section. 

The transform B'^{v) is defined implicitly in the standard manner using a 
Darboux transformation: one considers the following local gauge transformation 
of the Lax operator (lO) 


D^+i{u,v)Lj{u;tj;,tp*) = Lj{u;^,'tj}*)D+{u,v) , 


( 2 . 8 ) 


where the so-called Darboux matrices are assumed to be of the following 
form, 


D+{u,v) = 


V — uaj 
VCf 


-ubj 

—u 


(2.9) 


Here v is an additional variable related to the spectrality property of the Backlund 
transform yB+(w); see m and references therein for an explanation in the 
context of the Toda chain. N.B. the dressing matrix is singular for u = v, 
det Dj' {v, v) = 0, if we require that Uj = l — bjcj in close analogy with the case of 
the Toda chain discussed in loc. cit. The equality (12.81) together with the explicit 
form of the Darboux matrix (12.91) then fixes the map B^{v) : {ip, ip*) >—t {ip, ip ); 
compare with [28l Eqn (3.1)]. 


Proposition 2.2 The transformed variables {ipj,ipj) obey the following func¬ 
tional relations 

( Pp--ip.=v^j_^{l-ip*^j) 

{ ■ (2.10) 

In particular, the matrix elements in i2.9]} are Oj = 1 vip*ipj_i, bj = —vip*, 

Cj = 


Proof. A straightforward computation. Inserting (12.91) into (12.81) and compar¬ 
ing coefficients of powers in the spectral variable u yields the asserted equalities. 

■ 

The identities (12.101) allow one to compute the transformed variables {ip, ip ) 
via recurrence upon expanding them in power series with respect to the variable 
V. That is, if we set 


r>0 r>0 

then the above relations imply the recurrence identities 

r = +^j E 

a+b=r-l _ ( 2 . 12 ) 

= 1p*1p*+lipj,r-l 

Note that H+(0) = Id which fixes the initial conditions ip^ Q = ipj and ipj^Q = ip*. 
By construction the map B'^{v) leaves the Hamiltonians (12.Sp as well as (12.4p 
invariant, Tr{ip,ip*) = Tr{ip,ip ). Furthermore, one has [28l Prop 5]: 
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Proposition 2.3 B^{v) preserves the Poisson structure 

Therefore the transform B^{v) satisfies the defining conditions of an in- 
tegrable map |33| and it then follows from a general argument m that the 
transforms B^{x) and B^{y) commute, 

B+{x)oB+{y)=B+{y)oB+{x) . (2.13) 


One can also consider the inverse of the Darboux matrix (12.91) and one then 

.A. .^ 3(5 

arrives at a second transform B~(y) : ('0,^*) i—>■ ('0,'*/' ) which in case of the 
infinite chain is related to the inverse of B^{—v) defined via the relations [28l 
Prop 4, Eqns (3.13) and (3.14)] 


“ '^0 = 

^*3 - ^*3 = ^(1 - ^ 3 ^* 3)^*3 + l 


(2.14) 


The proofs are analogous to the previous case. Both systems of equations, (I2.10p 
and (12.141) . are approximations of the flow (12.71) . Similarly, one can investigate 
the transforms for the flow generated by —T„_i; see [281 Prop 2 and Prop 3]. 
These turn out to be related to the hermitian adjoints of the Q^-operators, we 
therefore restrict our discussion to (I2.10|) and (12.141) . 


3 Quantisation of the Ablowitz-Ladik system 

The following quantisation of the Ablowitz-Ladik model was first discussed by 
Kulish [T7] . 


dtPj — Pj+i — + Pj-i — P*jPj{Pj+i + Pj-i) 

dtp* = -P*+, + 2/3; - P*_, + P*PjiP*+3 + P]-i) 


where {Pj,P*j}jein now noncommutative variables satisfying the defining 
relation of the g-boson or g-Heisenberg algebra 'Hn{q)\ c.f. [13 Eqn (14)] and 
|3]. Throughout this article we assume q to be an indeterminate. 

Definition 3.1 Let LLniq) be the C{q)-algebra generated by {Pi, P*i,q^^'}*i=i 
subject to the relations 

r q^^P,=q-^^^p^q^^, q^^ P* = q^^^ P*q^^ 

{ P.P*-P*P.=S.,{l-q^)q^^^, P^P*-q^P*P^ = l-q^ ' 

N.B. the notation q^^ is purely formal, it does not mean that q^' is given 
by exponentiating another generator Ni although we will often write q°‘{q^')'^ 
as We have also made a small change to the usual definition of the 

g-boson algebra, see e.g. m, by multiplying one of the generators with an 
extra factor, f3j —>• (1 — q^)Py This allows us to derive from the second set of 
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relations in (13.21) the identity I3*jl3j = 1 — and, thus, we obtain the familiar 
quantisation formula [•,•] = + 0{h?) for the Poisson structure (I2.1|) 

when evaluating the indeterminate as g —>■ exp(i?i 7 / 2 ), 

[/3„/3;] = %(1 - g2)(i _ + 0{h^) , (3.3) 

compare with m Eqn (13)]. Similar to the classical case, /?* is initially an inde¬ 
pendent generator, but we may consider representations of the g-boson algebra 
where /3* is the hermitian adjoint of /3^-; see the Fock representation (13.141) be¬ 
low. This corresponds to the above mentioned reduction of the Ablowitz-Ladik 
chain to the NFS model and due to the change to imaginary time, t —> it, one 
then must choose ^ = ic with c G M being the coupling constant in the QNLS 
Hamiltonian. 

The form of the Hamiltonian closely resembles the classical one, 

n 

H - 2(1 - q^)N,), (3.4) 


where the iVj ’s are a set of additional generators not contained in the original 
g-boson algebra obeying 

iv,g^. =g^^iV„ /3,(iV,-%) = iV,A, P*{N,+5,,)=N,P* . (3.5) 


If we identify j3* with the creation and /3j with the annihilation of a g-boson at 
site i, then we easily recognise that Ni plays the role of a number operator on 
this site. We shall denote the total particle number operator by TV = 

In close analogy to the classical case (EH) one defines the monodromy matrix 
of the quantised system as a product of the local Lax matrices 


L{u) = Ln{u) ■ ■ ■ Li{u) 


A{u) B{u) \ 
C{u) D{u) J 


Lj{u) 


1 up* \ 
Pj u ) 


(3.6) 

and one then obtains the quantum versions of the spectral invariants EH of 
the classical system. 


T{u) = Tr L{u) = A{u) + z D{u) = ^ u^'Tr , (3.7) 

r=0 


where z is an additional quasi-periodicity parameter which we will need later 
on. The explicit form of the quantum integrals of motion has been derived in 
[T5l Prop 3.10, Eqn (3.43)]: set ai = PiP*i+i for i = 1,... ,n— 1 and a„ = zP^Pl 
(right or clockwise hopping) then 




(l-g2) 


r — 1 


(3.8) 


with [X, Y]q 2 = XY — q^YX and the sum is running over all cyclically ordered 
words w with letters 1 < ij < n occurring at most once. Moreover, Tg = 1, 
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Tn = z and we set Tr = 0 for all r > n and r < 0. Note that Ti = J2j=i 
and that Tn simplifies to T„ = Thus, as in the classical case one 

has a splitting of the quantum Hamiltonian (13.41) into ‘right movers’ and ‘left 
movers’. 

One now proves quantum integrability of the system (ETD via the following 
solution of the Yang-Baxter equation 


Ri2{u/v)Li3{u)L23{v) = L23{v)Li3{u) Ri2{u/v), (3.9) 


where the ii-matrix in End(C^ (8) C^) = EndC'* with the conventions used in 
the definition (|3.4I) of L reads 


R{u) 


/ 0 0 0 \ 

u—1 

0 0 

^ u—1 

0 ^ 1 0 

u—1 

\ 0 0 0 ^ ) 


(3.10) 


By the standard argument one obtains as an immediate consequence the follow¬ 
ing proposition. 


Proposition 3.2 The subalgebra An C Rn generated by the Tr is commutative, 

[Tr,Tr^]=0. 

Proof. One verifies that (13.101) is invertible. Using the cyclicity of the trace 
in (1X71) one obtains [T{u),T{v)] = 0. This proves the assertion for z = 1. The 
case of general z follows by the same argument using that R{u/v) commutes 
with (J °) O (J °). ■ 

We now consider a particular representation of the g-boson algebra which 
allows one to apply the quantum inverse scattering method or algebraic Bethe 
ansatz to construct a common eigenbasis. 


3.1 The Fock representation of q-bosons 

For the sake of completeness we recall the following representation of the g-boson 
algebra in terms of (/-difference operators, which can be found in Consider 
the ring of polynomials IH" = C((/)[^i,... ,^„] with rational coefficients in the 
indeterminate q, where the ^/s are some auxiliary variables. Let be the shift 
operator with respect to 

= /eSH" (3.11) 

and define Vi to be the following (/-derivative with respect to 

(P,/)(^) = _ (3,12) 

Si Q Si 
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In the limit <? —t 1 one recovers the ordinary partial derivative di. We introduce 
a hermitian bilinear form SH” x 91" —)> C((?) by setting 

(/,<?) = /(25i,-,25n)5(ei,-,Ulei=-=?„=o (3.13) 

and define = 91" to be the completion with respect to this inner product. 
We will simply write T for and we have that The following 

result is contained in m- 

Proposition 3.3 (Fock representation) The mapTLn — t EndJ^" defined via 

I3i^ 0-- (3.14) 

where is the multiplication operator with defines a simple module for Hn 
for all n > 1. 

There is an alternative formulation of the same module used in |151 Prop 
3.3]: let T C Hn be the left ideal generated by the elements /3j and (1 — 

Define a highest weight vector |0) = l+I, the “pseudo-vacuum”, in the quotient 
HnjT with 1 the identity element. For any partition A with at most n parts set 

l^) = nrlA- (3.15) 

(Q JmjiX) 

where mfiX) is the multiplicity of columns of height j in the Young diagram 
of A. It is easy to verify that the map |A) i-J- 0^=1 provides a 

module isomorphism Hn/T ^ and we shall henceforth identify both mod¬ 
ules. We denote by (A| the dual basis of (13.151) and the map .F" —>■ F"" given 
by (^1 b\{q^)\X) with F" denoting the dual space and b\ = Y\j>i{q^)mj{x)^ 

then introduces an inner product on F". Note that if q is evaluated at a root 
of unity the module ceases to be simple. 

Remark 3.4 The eigenspaces of the particle number operator N = 
are the subspaces Tff C F" spanned by the vectors {|A)}Aj=fc for k € Z>o. In 
m it has been shown that F^ carries an Uq{sln)-action and can be identified 
with the Kirillov-Reshetikhin module of highest weight kuii with uii being the first 
fundamental weight. The basis iS.15\) is then Lusztig’s canonical or Kashiwara ’s 
global crystal basis for this module. 

3.2 Completeness of the Bethe ansatz 

The existence of the highest weight representation (13.141) allows for the appli¬ 
cation of the quantum inverse scattering method [5] or algebraic Bethe ansatz 
and one obtains the following important result. 

Denote by k = C{{g}} the algebraically closed field of Puiseux series in the 
indeterminate q. 
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Theorem 3.5 In the Fock representation the difference operators given by hS.Sil 
possess a common eigenbasis {|2/a)}a where 

(3.16) 

and A ranges over the partitions with at most n parts. The partition A labels the 
solutions y\ = {yi,..., y\^) G of the following coupled set of equations 


yf = z 


n 




- Vj 

yi - yjo^' 


z = 1,2,..., Ai . 


The eigenvalue equation reads, 


T{u)\yx) 


f 1 - uq^yj 

\i=i 1 ~ 




I^a) ■ 


(3.17) 


(3.18) 


Remark 3.6 The quantum inverse scattering method for the q-boson model 
was discussed in and m where the Bethe ansatz equations can 

be found. A discussion of the coordinate Bethe ansatz for the discrete QNLS 
model arriving at the same set of equations appeared also in JSOf . where a proof 
of the completeness of the Bethe ansatz solutions can be found for q = e with 
— 1 < e < 1. Completeness of the Bethe ansatz for q an indeterminate together 
with a coordinate ring description of their solutions as quotient of the spherical 
Hecke algebra was obtained in Section 7j. In order to find solutions of Ili3.17\ ) 
one then needs to work over the algebraically closed field of Puiseux series in q. 


4 Two Q-operators for the g-boson model 


In order to construct the quantum analogue of the Backlund transform (I2.10p 
and its inverse one considers solutions D^, to the Yang-Baxter equation 

Di2iu,v)Li3{u)Lf^iv) = Lf^{v)Li3{u)Df2iu,v), (4.1) 

where L{u) is the local quantum Lax matrix in (13.61) and D^{u,v), L^{v) are 
respectively 2 x ^ and ^ x ^ matrices with entries in the g-boson algebra. 
As explained in for the Toda chain, the Yang-Baxter equation (14.11) should 
be seen as a quantum analogue of the classical relation (j2.8p for the Darboux 
matrix: due to the noncommutative matrix elements in the quantum case an 
extra matrix is required. Thus, in the quantum system one faces the problem 
of hnding two matrices, and L^, to construct the quantum analogue of the 
Backlund transform (12.101) . 

Define the following half-infinite matrices with entries in the g-boson algebra 
Hn (c.f. [m Eqn (3.20)]), 


Lit (^') = i-^y 


MrPT 




(4.2) 


12 














and 


(4.3) 


L-{v) = 


m(m+l) PT {PIT 




m,m'>0 


The operators which are not present in the classical equation enter in the 
definition of the Q^-operators 


Q^{v) =TT:z^L^{v)---L^{v)Lf{v) 


(4.4) 


r>0 


where = (z^ ) and the trace is formally defined with Q^(v) G C|u](8)'H„. 

That is, Q^(v) should be thought of as current operators: as formal power series 
in the variable v with the coefficients being elements in the q-boson algebra 

Un- 


Lemma 4.1 The coefficients in the expansions U-4]) are 

(/3D“"(/3i^;)“^---(/3„i^:)“"-/3“ 


g+ = (- 1 )’'^ 

ahr 




(4.5) 




cuhr 


{q^U---{q^)o.^.Aq^) 




i=l 


where the sums range over all compositions a = (ai,..., a„) G Z>q of r > 0. 

Proof. Via induction in n one verifies that the monodromy matrices T^{v) = 
L^iv) ■ ■ ■ L 2 iv)Li(v) are given by (c.f. [HI Eqn (3.48)]) 






iqAr 


(9^)01•(g^)a„_i 


and 


^m'miy) = 

„ iq^)m (<?^)ai---(g^)a„_i(g^)a„ 

where the sums range over all compositions a = (oi,..., On-i) G Z"q^. Sum¬ 
ming over the diagonal matrix elements with m = m' and fixing the degree of 
V the assertion now follows. ■ 

4.1 Functional relations in the g-boson algebra 

The characteristic property which prompts us to identify (14.41) as the analogue 
of Baxter’s Q-operator for the g-boson model is the following set of functional 
relations. The first was already proved in [151 Prop 3.12], we repeat it here for 
completeness, the second identity (14.81) is new. 
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Theorem 4.2 We have the following identities in the q-boson algebra 

T{u)Q^{u) = Q^{uq^) +/\{u)Q^{uq~‘^), (4.7) 

Q~{u)T{u) = Q~[uq~'^) A-/^{uq^)Q~ {uq^) , (4.8) 

where the coefficient is given by A(u) = zq^^u^. 

Proof. For the proof of the first equation we refer the reader to nn Prop 
3.12]. The proof of the second equation follows along similar lines. Namely, we 
consider the action of Lf^{u)L 23 {u) on the tensor product V (Si Hi with 
^ = ©™>o Cvm, where we identify the basis vectors in as vo,vi. Let X be 
any element in Hi. Then a straightforward computation yields 


L^^{u)L23{u)Vm ®Vo® X = 


ymgm(m+l) 


m'>Q 


as Vo® r\ii*rx+ vm' ®vi® r'{p*ri3x 


and for m > 0, 


Li^{u)L23{u)Vm-l ®vi® X = 

m m{m-l) , , 

- 5] Vrr,! ®Vi®r + Vm' ® Vo ® T 

W jm-1 

Consider now the subspace W C F C) spanned by the vectors {wm = Vm® 
+Wm-i I m > 0} and wo = vo®vo. Employing the commutation relation 

= (/3*)™-i(l - q-2™) + g-2™/3(/3*)™, 


which is easily verified by induction, one then finds from the two identities above 
that on W ® Hi we have 

„,m -2m m(m+l) _ , 

Lf^iu)L 23 iu)Wm ®x= - E ® • 

(q jm-1 m'>0 

Next we consider the action on the quotient space W' = F C> jW. A basis is 
given by the vectors {vm ®vi | m > 1} and one arrives at 


L^^(u)L23{u)Vm ®vi® X = 


ym+l qm(m+l) 


{<P)r. 




m '>0 

where the omitted terms lie in IF 0 'Hi. Since 


m ^m'-\-l m+1 _ ^2m+2 m^2N 
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we deduce that the product of L-operators block-decomposes as follows, 
ir3(«)^23(u) = ^ ^ ) 

and the assertion (14.81) is now easily obtained by taking the trace of the product 


Tq {u)To^{u) = Lo „( u ) Lo ' n ( M ) • ■ ■ Lq^{u)Lo^2{u)Lq^{u)Lo,i{u) 


of the monodromy matrices and applying the above block decomposition to each 
factor LQj{u)Lo>j{u). ■ 

In addition to the functional equations Kn and (|4.8I1 one needs to prove 
that 

[T(u), Q±(u)] = [Q±(u), Q±(u)] = [Q±(u), QT(u)] = 0 (4.9) 

for arbitrary values of the variables u,v. The first and second relation for the 
Q+ operator has been proved in m Props 3.7 and 3.8] by constructing explicit 
solutions to the Yang-Baxter equation (HU) with [m Eqn (3.23)] 


D'^{u, v) 


1 _ HqSAT \ 

V I 

/3 —u/v j 


(4.10) 


as well as finding an additional solution of the equation i')^i 3 (w)^^ 3 ('*^) = 

(v)L^^(u)Ri 2 {u,v). The above expression for D~^ yields the desired 2 x ^ 
matrix in if the g-boson algebra elements in are evaluated in the 

Fock space representation (13.1411 . Note the close resemblance of the quantum 
Darboux matrix (14.1011 with the classical one (12.911 . In principle one can proceed 
in the same manner for Q~ and one then finds 

D-{U,V) = ^ _ fj27V+2 ) (4.11) 

which gives a second quantum Darboux matrix which closely resembles the in¬ 
verse of ()4.10ll [m Eqn (3.25)] in accordance with the classical case. To estab¬ 
lish (14.911 one needs to find yet two other solutions of the Yang-Baxter equations 
Ri 2 iu,v)L^ 3 {u)L 23 {v) = i^ 3 (w)LC 3 (u)i?C 2 (u,u) and R[ 2 {u,v)L+^{u)L 23 {v) = 
L 2 ^{v)L^^{u)R'i 2 {u,v). Here we shall instead take advantage of the already es¬ 
tablished functional relations (14.71) . (14. 8p as well as Prop [XU to give a much 
shorter and less computational proof of (14.91) . 

Corollary 4.3 The operator coejficients of the formal power series T{u), Q^{u) 
all commute, i.e. we have \Tr, Qf] = [Qt^ Qj] — ^ r, s > 0. 


Proof. Rewriting the functional equations (14.71) and (14.81) in terms of coefficients 
we find that 

r 

(1 - , 

S^l 

(q-2’--l)Q- = + (4.12) 

S = 1 
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where we set = 0 for r < 0. Noting that Qq = 1 by definition and that for 
each r > 0 the terms on the right hand side of both identities only involve Qf 
with s < r we find that the Qf are polynomials in the TQs. But since the latter 
commute among themselves, [Tj.,Ts] = 0, the assertions now easily follow. ■ 


Corollary 4.4 The transfer matrix T{u) and Q^{u) obey the additional rela¬ 
tions 


l = Q^{u)Q {uq - zvTq^^Q^{uq '^)Q (u) (4-13) 


and 


T{u) = 


Q+{uq^) A{u)Q~''{uq 

A{uq^)Q~ {uq^) Q~{uq~'^) 


(4.14) 


As explained above these equalities should be understood as identities in terms 
0fTr,QfGHn- 


Proof. Consider the triple product T{u)Q'^{u)Q (u). From (|4.7I) and (14.81) one 
deduces that the formal power series W{u) = X]r>o = Q'^{u)Q~[uq~'^) — 
Q+ {uq~‘^)Q~{u) obeys W(u) = W{uq^) and, hence, Wr = 0 for r > 0 
since q is an arbitrary indeterminate. Setting u = 0 we hnd bF(0) = 1 and the 
first assertion follows. The second equality (14.141) is now easily deduced from 
(|4.13|) by using once more (14.71) and (14.81) . ■ 


4.2 Specialisation to the Fock space representation 

It is worth emphasising that so far we have worked on the level of the q-boson 
algebra, that is, the functional relations derived hold regardless of the represen¬ 
tation chosen. We now specialise to the Fock space representation (13.141) . 

In the Fock representation the coefficients (Ii3]) . (ITO of the Q-operators 
become the following g-difference operators. 


Qf 

Qf 


ahr 


(-1)^(1-q^yJ2 


'' ck/ 


f n 

-n-'° 


ahr 




(4.15) 

(4.16) 


where ai = (a„, oi, 02 , • ■ ■, ctn-i)- We can now take matrix elements of the 
Q^-operators and the L^-operators in the Fock space representation. For the 
local Lax matrix (13.61) we obtain for the vertex configurations in Figure [T] the 
Boltzmann weights [m Fig 3] 

(0, m|L(M)|0, m) = (1, m — 1 |L(m)| 0, m) = 1 (4-17) 

and 

(0, m -I- 1|L(m)| 1, m) = m(1 — (1, m|L(u)|l, m) = u . (4-18) 
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From the L^-operators we find the following Boltzmann weights for the vertex 
configurations in Figure [2l 


{c,d\L'^{u)\a,h) = (— 


(4.19) 


and 


(c,d|L-(u)|a,6) 



(4.20) 


with It -for m > n and zero otherwise. The matrix elements of 

the (5^-operators then yield the weighted sums over the vertex configurations of 
a single lattice row, the row-to-row transfer matrices. According to (|4.9II these 
define exactly solvable lattice models. As we have interpreted Q^{u) so far as 
formal power series with coefficients in the g-boson algebra we briefly explain 
how these give rise to well-defined operators in Fock space. 

For fixed particle number k the formal power series (5+(u)|^n is a well- 
defined operator as = 0 for r > fc. In fact, from (14.511 one deduces that 

Q+Ijrn is the discrete translation operator: the sum in (14.511 ioi r = k is the 
sum over all possible fc-particle configurations shifting each particle by one site 
forward. If r > fc the corresponding operator Q+\jrn would shift more particles 
forward than are in the system, hence its matrix elements vanish. 

In contrast, the coefficients of Q~iu)\j^n are in general all nonzero on the 
fc-particle space, since according to (j4.3ll and (j4.8ll particles are first created at 
the neighbouring site before they are annihilated at their place of origin. Thus, 
under the action of Q~ iu)\j^n a single particle can travel several times around the 
lattice picking up a factor z each time it completes a round-trip. We therefore 
interpret the power of the quasi-periodicity variable z in (13.711 as a winding 
number. For any finite winding number p > 0 the corresponding coefficient of 
z^ when expanding Q~{u)\jrn as a power series in z, is a well-defined operator, 
since we have now limited the number of round trips a particle can make and 
there are only a finite number of particles in the system. 


Corollary 4.5 In the Fock space representation we have that the Bethe vectors 
id. 1 (A) form a common eigenbasis of {Tr,Qf} with 

Q^iu)\y\) = ri (1 - M Vj) \yx) . (4.21) 

1=1 

In particular, the eigenvalues of Qf are the elementary symmetric functions in 
the Bethe roots yi. The eigenvalues of Qf are then derived via {4.13\) . 

Proof. From (14.121) it follows that the Bethe vectors (13.161) are eigenvectors of 
as the latter are polynomial in the T^’s. The identity (14.71) together with 
(j3.18|) gives the eigenvalue Q^{u,y\) — njii(l ~ Vj) in <14.211) . ■ 

Naively one might expect that Q~{u) has analytic eigenvalues as well. Let 
us evaluate u in a neighbourhood U C C of m = 0 and for a fixed number of 
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lattice sites n and particle number A: = Ai set z in (13.71) to z = e{n,Xi) with 
0 < £ <C 1. If one can show that the following expression converges, 


F{u,yx) = Q+(m,2/a)^ 

r>0 


^rn ^nr{r-\-l) ^2rXi 

Q+ {uq^^ ,y\)Q+ (uq' 2»'+2 ^ ) 


(4.22) 


then the Bethe ansatz equations (13.171) imply that the residues at u = q~'^"^/yj 
vanish and, hence, that F is analytic in u. Via analytic continuation in e one 
then defines F outwith the region of convergence. The identity (j4.8l) would 
then imply that F{u, y\) = Q~(u, y\) is the eigenvalue of Q~{u). However, the 
expression (14.221) is also a power series in q with the Bethe roots being potentially 
Puiseux series in q. In order to control the convergence of the expression (I4.22|) 
one needs to know the dependence of the Bethe roots on q which is a difficult 
and technical issue. We hope to address this problem by different means, which 
go beyond the discussion in this article, in future work. 

Remark 4.6 After the construction of Q'^ in m Section S] an alternative 
expression for a Q-operator has been put forward m in the Fock representation. 
In loc. cit. the Q-operator is .stated in terms of a kernel function for the Jackson 
integral (Eqns (54) and (62)) which requires special convergence conditions, 0 < 
q < 1, and an upper bound on the growth of functions in the state space 
The derivation of this kernel function in loc. cit. postulates the existence of a 
q-analogue of the S-function for the Jackson integral and I was unable to verify 
whether the constructed operator is either of the two operators constructed here. 

N.B. in the proof of (14.211) we have reversed the usual logic and derived the 
spectrum of the Q'^'-operator from the spectrum (13.181) of the transfer matrix 
and the functional relation (HTTl) . As the system is solvable via the Bethe ansatz 
and completeness has been proved [m Section 7], the significance of the Q~^- 
operator in the present context is not that of a mere ‘auxiliary matrix’ as in 
Baxter’s original work [5] (and references therein) where it is used to find the 
physically relevant eigenvalues of T. Instead the Q'''-operator acquires in our 
setting a direct physical significance as it describes the discretised time evolution 
(IQ) of the quantum Ablowitz-Ladik chain, as we will see next. 


5 Quantum Backlund transformation 

We define the quantum analogue of the Backlund transformation in the 

Fock space representation by setting B'^{v) : End J"” —>■ C[u] 0 End J"" with 
O 1-4 {v)OQ^{v)~^. The quantum Backlund map is well-defined because 
of the following result. 

Proposition 5.1 In the Fock representation the inverse of the Q'^-operator 
exists and is given by 

Q+(u)-i = ^?;’'det((-l)i-*+^Q)L,+j)i<*j<r , (5.1) 

r>0 
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where the coefficients for r < n simplify according to 

, r<n. (5.2) 

Proof. The first identity for the inverse operator follows from the Bethe ansatz 
result: in the eigenbasis (13.161) the inverse of Q'^{v) is given by the diagonal 
matrix with entries 

j=l ^ r>Q 

with hj. denoting the complete symmetric functions [B Ch. I]. Using the known 
determinant relation between elementary and complete symmetric functions, 
hr — dei{ei-i+j)i<i^j<r, the first assertion follows from (14.211) . The second 
assertion is now an immediate consequence of (14.1311 when expanding in the 
variable v. ■ 

We are particularly interested in the image of the g-boson algebra generators 
under the quantum Backlund transform, 

Pj{v) = Q^{v)l5jQ^{v)~'^ = ^v'^Pjr ■ (5-3) 

r>0 

We define Pj{v) and fijr in nn analogous fashion. N.B. j3j{v) {Pj{v))* since 
Q+ is not unitary in the Fock space representation. For ease of notation we will 
often suppress the explicit dependence on the variable v but the reader should 
keep in mind that I3j{v) and j3j{v) are power series in v with coefficients in 
EndJ^" as we are now working in the Fock space representation. Clearly, the 
transformed quantum variables still obey the g-boson algebra relations 

and, thus, we obtain a one-variable family of representations of T-Ln- Moreover, 
by construction and because of (14.9p the quantum integrals of motion are left 
invariant, 

TrCHj.Kj) = Q+{v)Tr{P^,P*)Q+{v)-^ = Tr{l3^,fi*) . (5.4) 

5.1 The image of the quantum Backlund map 

We now derive a set of functional relations for the transformed quantum vari- 
ables and show that they are an exact match of the classical relations 

(|2.10l) which describe the discretised time flow (12.71) . 

Theorem 5.2 The quantum Backlund transformed variables {/3j{v), fij (z;)} obey 
the functional relations 

r - (3^ = n(l - 

1 ■ (5.5) 

These relations determine {fij{v),fij{v)} via recurrence when expanding in the 
variable v. 
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Proof. A somewhat lengthy but straightforward computation. We give the 
intermediate steps and leave some details to the reader to verify. Set aj = 
/3j/3j+i- Via induction in m one establishes the identities 

[a-i, /?,-] = - (1 - , 

From the explicit expression (14.51) one then derives the recurrence relations 

[Q+, ( 3 *] = -/3;+iQt - /3;+i ■ 

Multiplying with u’' and subsequently summing over r on both sides of these 
two equalities produces 

[Q+(u),/3^-] = 

[Q+(u),/3;] = -vp*^,Q+{v)q^^^ -vp*^,[Q+{v),/3*]f3^ . 

After multiplying with from the right we obtain (recall that q'^^^ = 

1-W 


r = u(i - P*P,)P,., - vp*[p, - 

1 ~p* - p* = vp*^,Cp*~p^ -1) - vp*+,{p* - p*)~p^ 

from which the desired equalities in (15.51) now easily follow. ■ 

We demonstrate that (15.511 allows one to compute the image of the g-boson 
generators via recurrence. Namely, making an analogous power series expansion 
Pj(y) = J2r>o'^^^j,r the classical case (12.lip we find for the first few 

coefficients Pj q = Pj and 

^,,3 = /3,_3(1-/3;-2/3,-2)(1-/3;-i/3,-i)(1-/3;/3,) 

-2/3,_2/3,_i^;(1 - - P*P,) ■ 

Similarly, we find for Pj that ,5^ g = P* and 


Pi, = -(1-P*P,)P*^, 

Pi, = P^_,P*P*^,il-P*P^) 

Pi, = -pUp- ^/3;+i(i - P*P,)+ /3,_2/3;/3;+i(i - pup,-,){i - p*p,) ■ 

We can interpret the above formulae as describing discrete time steps under the 
evolution (12.7p : since the quantum relations (|5.5p and classical relations (12.101) 
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~ ~ 5)t 

are the same, these expressions yield the corresponding coefficients and 

~ ~ ~ 5(! ~ jf: 

in (I2.12P upon replacing /3^,, —)• and j3j j. —>■ 4’j,r- 

Naturally, one wants to extend the discussion to include the (5“-operator 
and one then needs to show that its inverse exists. However, when proceeding 
along the same lines as for the (3“''-operator one must first prove convergence of 
the supposed eigenvalues (14.221) . Since the latter is currently an open question, 
we state the following preliminary result. 

Lemma 5.3 We have the following commutation relations between the genera¬ 
tors of the q-boson algebra and Q~[u) 

f Q~iv)l3j - l3jQ-{v) = -vl3j_i [(5“(u) +/3j(3“(u)/3*] 

\ Q-iv)p*-p*Q-iv)=v[Q-iv)-p^Q-iv)/3*]p*+, 

Proof. Via a similar computation as in the case of Q'^ one shows the commu¬ 
tation relations 


[Q-,/3^] = , 

after multiplying with v'^ and summing over r the assertion follows. ■ 

Define a one-variable family of operators {fij{v),fij{v)} via the functional 
relations (compare with (12.141) 1 

As in the previous case (1531) the relations (E3 also determine {^j(v), (v)} 

via recurrence when making the analogous power series expansions /fjiv) = 
Pjw' Under the assumption that Q~{v)~^ exists, it then follows from 
(15.61) that fij{v) = Q~{v)fijQ~{v)~^ and ^j{v) = Q~{v)fi*Q~{v)~^. 


6 Prom Backlund transformations to 2D TQFT 

We now link the quantum Backlund transform (15.51) to the 2D TQFT con¬ 
structed in [151 Section 7]. Consider multivariate deformations of the g-boson 
algebra by setting 

...,xe)-.= Pj (6.1) 

with i < n. We define I3j{xi,...,xe) analogously. We shall concentrate on the 
Q’^-operator because both solutions are related via the functional relation 
()4.13p . which shows that the transformation induced by Q~ can be constructed 
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from (5+. Since the Q^{xiys in (ICT commute with each other - reflecting 
the analogous property (12.131) of the classical Backlund transform - one can 
expand the product of Q^-operators In any basis In the ring of symmetric func¬ 
tions in the variables Xi. Set P\{xi,... ,xi;q) = P\{xi,... ,xi;q'^,0) where 
P\{xi,... ,xe;q,t) are the celebrated Macdonald functions [T^. Then the ex¬ 
pansion 

n Q^{xi) = ...,xe;q) (6.2) 

i=i ^ 

with the sum running over all partitions with at most I parts defines uniquely a 
set of commuting polynomials {Qa} in the g-boson algebra "Hn ® C| 2 ], where z 
is the “quasi-periodicity parameter” in (13.71) which we treat as formal variable. 
Note that the operators Q\ can be defined explicitly as polynomials in the Q^’s, 
see [m Def 3.3]. For the discussion in this article their Implicit definition via 
(16.21) sufHcefl 

In the Fock representation it follows from the results in [151 Section 7] that if 
A has a column of height n then Q\ = where A Is the partition obtained 

from A by removing all columns of height n and mn{X) is the multiplicity of these 
columns. Setting z = 1 (periodic boundary conditions) we therefore restrict 
ourselves to i = n — 1. We recall the following results from (THl Thm 7.2, Lem 
7.7 and Cor 7.3 ]. 

Theorem 6.1 (1) The {Q 3 ,} form a basis in the ring An of quantum integrals 
of motion generated by the Tr’s. (2) Let A, /i, he partitions with at most n parts 
and fii = vi = Xi = k G Z>o. Then 

Q\Qh = = (^IQaIm) , (6-3) 

I'l—k 

where the expansion coefficients are the fusion coefficient of a 2D TQFT with 
A^~'_(0) being the SU{n) WZNW fusion coefficient at level k. If pL^ or 

Pi Ai the matrix element is zero. 

6.1 A brief summary of 2D TQFT 

For the sake of completeness and to make this article accessible to a wider audi¬ 
ence we briefly summarise the definition of a 2D TQFT, for the precise definition 
we refer the interested reader to the abundant literature on the subject; see e.g. 
the text book [la¬ 
in modern mathematical language a 2D TQFT, Is a symmetric monoldal 
functor Z : (2Cob,U) —>• (Vectk,®) from the category of 2-cobordisms Into 
the category of finlte-dlmenslonal vector spaces over some base field k. The 
objects in the category 2Cob are circles closed strings, and we define a 

^The reader should note that the definition of the Q+-operator in this article corresponds 
to the operator G'{—u) in 1151 Prop 3.11, Eqn (3.46)] which explains the extra sign factor in 
(ES compared to Eqn (3.56) in loc. cit. 
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Figure 4: Depiction of elementary 2-cobordisms. From left to right: the identity 
map Id : the ‘pair of pants’ m : Z(S^) 0 ^(S^) —>■ ^(S^), 

the ‘elbow’ : 2’(§^) 0 Z’(S^) —>■ k, the ‘cup’ 1 : k —Z(S^) and the ‘cap’ 
tr : Z(§i) ^ k. 


product on them via the disjoint union U. The functor Z maps each closed 
string onto a vector space Z(S^) while preserving the product structure, i.e. 

U • • • U §^) = Z’(S^) ig) • • • 0 Z’(§^) where ig denotes the ordinary tensor 
product of vector spaces. One includes the case where the circle shrinks to a 
point pt and sets .^(pt) = k. This explains the monoidal property of Z. Asking 
Z to be symmetric refers to the fact that each circle carries an orientation: let 
be anti-clockwise oriented and denote by the clockwise oriented circle. 
Then one demands that a change in orientation in 2Cob corresponds to taking 
the dual space in Vectk, i.e. Z(S^) = Z(S^)*. 

The morphisms between objects in 2Cob are two-dimensional compact man¬ 
ifolds which interpolate between two sets of circles, one representing the “in¬ 
states” U • • • U and the other the “out-states” U • • • U in physics ter- 

m n 

minology. The morphisms in 2Cob are called 2-cobordisms and can be thought 
of as a 2-dimensional analogue of Feynman diagrams with m incoming particles 
and n outgoing ones. Each such 2-cobordism is mapped under Z to an element 
in Hom(Z(S^)®'", i.e. a linear map between two tensor products of 

the vector space Z(§^). 

Because of the functorial property Z is fixed by specifying its values on 
certain “elementary” 2-cobordisms shown in Figured which introduce on Z{S^) 
a product m : Z(S^) (gi Z(§^) —>■ Z(S^), an identity element 1 : k —>• Z(S^) and 
an invariant bilinear form (-j-) : Z(S^) g Z(§^) —>■ k. Instead of the latter, one 
often considers the trace functional given by tr(-) = (1|-) which defines a map 
—>■ k. These maps endow with the structure of a Frohenius algebra 

and the latter is called symmetric if the product is commutative. 

6.2 Prom Q-operators to TQFT fusion matrices 

Hence, in order to prove that (lOl) defines a 2D TQFT we need to endow 
the algebra An generated from the quantum integrals of motion (13.81) with the 
structure of a symmetric Frobenius algebra. From (|3.8p and (j4.5ll . (14.61) we infer 
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that An block decomposes over the subspaces of fixed particle number, 

An — An,k 7 An,k C End *7^^ 

k>0 

with dim = dim^„_/c = Fix fc G N and set k = Cflg}}, the field of 

Puiseux series in q. Consider the vector space ^fe(S^) = An,k together with the 
matrix product (16.3p in An,k C End and 

lfc:=Q 0 and , (6.4) 

where A = {k - Xn-i, ..., /c - A 2 , A: - Ai) and 63 ^ = 

being the multiplicity of columns of height j in the reduced partition A. 

Theorem 6.2 ( |15L Thm 7.2 & Cor 7.3]) The maps \6.4\) together with li6.S\} 
define a commutative Frohenius algebra, that is Z},, : (2Cob, U) —>■ (Vectk, ®) with 
^fe(S^) = An,k is a symmetric monoidal functor. 

Remark 6.3 The fusion coefficients i6.3\} can be explicitly computed using the 
algorithm in m Section 7.3.1]. The latter is based on an algebra isomorphism 
which maps An,k on a quotient of the spherical Hecke algebra; see f73l Thm 7.3] 
for details. Alternatively, one can extract Af^^{q) from the combinatorial defi¬ 
nition of cylindric generalisations of skew q- Whittaker functions: taking matrix 
elements in hd.A) one obtains symmetric functions Pxjd/fi via 

_^ n—1 

^z‘^P\/d/^t{xi,...,Xn-i;q) = {X\ Q^{xi)\ff) ■ 

d>0 

Here Xfdfp, denotes a cylindric shape, a periodically continued skew diagram, 
and each Px/dj^i equals a sum over q-weighted cylindric semi-standard tableaux; 
see m Section 6.2]. 


6.3 The quantum Backlund transform in terms of fusion 

The following result connects the fusion coefficients (16.31) for the simplest, one- 
variable case, with the quantum Backlund transform (15.5|) . 


Proposition 6.4 The solution of the quantum Backlund transform id.51) satisfy 
the following linear systems involving the fusion coefficients nnt i. 


a-\-b—r p 


E E 

a-\-b—r p 


1 _ q2mj(fj,)+2 







(9) , 


(6.5) 

( 6 . 6 ) 


where the notation /3jpL and ]d*p, stand for the partitions obtained by respectively 
deleting and inserting a column of height j in the Young diagram of p,. 
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Proof. This is an immediate consequence of the definition (16.21) of the fusion 
matrices and (1631) which yields On fti® other hand the 

definition (15.31) implies $j{v)Q'^{v) = Q'^{v)f3j. Taking matrix elements in the 
latter equality and then expanding in the variable v yields (16.51) . The proof of 
(1631) follows along the same lines observing that /3j|^) = (1 — 
according to (13.151) . ■ 

For small particle numbers k there is a third way to compute the fusion 
coefficients (j6.3p using the following “operator-state correspondence” in the Fock 
space [151 Cor 7.3]: let Ifc") = |0) denote the basis vector (I3.15P with k 

particles sitting at site n and all other sites being empty. Then the basis vectors 
(|3.15p with Xi = k are obtained by |A) = where A is the partition 

obtained from A by removing all columns of height n. Thus, we can introduce a 
fusion product on the fc-particle space by setting |A) * |/r) = QxlfJ.) and the 
resulting algebra is isomorphic to the TQFT (16.3p . In order to obtain the fusion 
coefficient we then need to explicitly compute as a polynomial in the Q+’s 
and let them act on the state vector j/i). In the final step one then removes 
all columns of height n in the partitions in the expansion of as these 

correspond to the trivial representation of SU{n). We demonstrate this on a 
small example. 

Example 6.5 Set n = 3. We wish to relate the fusion coefficients i6.3\) for 
particle numbers fc = 2, 3 via the relation In the latter equation choose 

j = 1 and pL = (2, 2,1) = ^ . Thus, j3\pL = and the reduced partition 

with columns of height n = 3 deleted is simply /i = 0 . The state \pt) is a 
2-particle state with one particle sitting at site 2 and the other at site 3. The 
state flip) = (1 — q^)|3, 2,1) is a 3-particle state, with one particle sitting at 
each site. According to 16'. 61) we have the identity 

(1 - = (A|^;, 2 |m) - (6-7) 

P P 

with fi Q = f*, fii = o,nd P 12 = T’/ie first term on the 

right hand side in the above equation gives 

Km = = (1 - 9")(i - K\^) ■ 

That is, the first matrix element in is nonzero only for A = (3, 2). For the 
other terms we need to compute the 2-particle fusion coefficients and 

first. This can be done using the operator-state correspondence |A) = Q 3 ^| 3 , 3) 
or the algorithm in m- One finds (we simply write the Young diagrams of the 
partitions instead of the ’s) 

m * 3 = □ and □ * g = P + (1 -f 0 . 

Setting q = 0 one can verify that these give the correct WZNWfusion coefficients 
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of SU{3) at level k = 2. Thus, we obtain 


p 


-(l-Q^)(l + g2)| 


)-(l-9Vlff) 


( 1 - 9 ")!^) 


for the second and third term on the right hand side of Adding all three 

terms we obtain the 3-particle fusion coefficients i) thus, arrive at the 

3-particle fusion product (after deleting columns of height n = 3 in all partitions) 

m * P = (1 + (□ + B + 

Again one can check that in the limit q — 0 these are the correct SU{3)-WZNW 
fusion coefficients at level k = 3. As demonstrated on this simple example one 
can use the linear system to compute fusion coefficients recursively. 


7 Conclusions 

In this article we have added new insight to the quantisation of the Ablowitz- 
Ladik chain (HID by using Baxter’s concept of the Q-operator. 

7.1 Summary of results 

The classical system (HID possesses a set of integral of motions, a Poisson com¬ 
mutative subalgebra , generated by the trace of the Lax operator (12.31) . 
{Tr,Tri} = 0. These integrals of motion generate additional “elementary” flows 
which can be used to construct the physical time flow given by the Hamilto¬ 
nian (j2.2ll . Suris constructed in [28] discrete approximations of these elementary 
flows in terms of Backlund transformations given in terms of a set of functional 
relations (|2.10l) , (|2.14|) . In [17] the q-deformed oscillator or q-boson algebra 'Hn 
was considered as a quantisation of the Poisson algebra (12.11) underlying the 
classical chain. 

In this article we have discussed the quantum analogues of the classical 
Backlund transformations using the g-boson algebra. We have stated the explicit 
construction of two infinite families {Q^jr-eN of elements in the g-boson algebra, 
see dUD and (|TO . which commute among themselves and with the quantum 
integrals of motion An generated from the higher Hamiltonians of the 

quantised Ablowitz-Ladik chain. Introducing the associated current operators 
Q^(y) we derived their algebraic dependence in terms of a set of functional equa¬ 
tions: Baxter’s TQ-equation (|4.7p and (14.8p . which is a second order difference 
equation, and the related Wronskian (14.131) showing that is closely related 
to the inverse of in the limit of large lattice sites. We showed that these rela¬ 
tions hold on the level of the g-boson algebra without the need to specify a rep¬ 
resentation. Once we fixed the Fock space representation (13.141) we were able to 
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show in addition that Q'^{v) is invertible for the periodic chain by employing the 
completeness of the Bethe ansatz for the q-boson model proved independently 
in [50] for — 1 < q < 1 and in m for q generic. This allowed us to consider the 
one-variable family of similarity transformations O i—?> Q'^{v)OQ'^{v)~^. Ap¬ 
plying these transformations to the generators of the (noncommutative) q-boson 
algebra (15:51) we have shown that their image obey a set of functional relations 
(1531) which match exactly the relations (|2.10l) of the classical Backlund trans¬ 
formation. This is an unexpected result: the discrete time flow of the 

classical and the quantum variables is the same. 

An open question is the extension of this result to the Q~ (u)-operator for 
which further investigations are needed to establish the existence of its inverse. 
As a preliminary result we showed that assuming the existence of Q~{v)~^ a 
second set of functional relations (EZj) follows which again match the classical 
relations (12.141) associated with the inverse of the classical Backlund transfor¬ 
mation. 

In the final section we then considered multivariate versions of the Backlund 
transformation via composition. The novel result is that they generate the 
fusion matrices of a 2D TQFT constructed in [15]. This TQFT is a q-deformed 
version of the SU{n)k-'WZNW fusion ring or Verlinde algebra which is recovered 
in the strong coupling limit q —^ 0; the related combinatorial description of the 
WZNW fusion ring and its relation to quantum cohomology can be found in 
m and m- Because of this limit, it was suggested in m Section 8] that this 
TQFT might be related to the one considered by Teleman [29] in the context of 
twisted AT-theory. 

7.2 Integrability and topological quantum field theory 

We explain how our findings differ from previously known connections between 
classical and quantum integrable systems and TQFTs starting with the quantum 
case. 

Nekrasov and Shatashvili developed in a series of works the following con¬ 
nection between TQFTs and the Bethe ansatz equations of quantum integrable 
models m- starting from a 4D Yang-Mills theory one can consider its vacua 
in the infrared limit. The moduli space of vacua is described by a TQFT in 
terms of a set of equations matching the Bethe ansatz equations of a quantum 
integrable model using the Yang-Yang functional. For the q-boson model this 
approach has been followed recently by string theorists in [35] and [10] who put 
forward generalised versions of WZNW and Chern-Simons field theories. 

In contrast our construction of a TQFT using Baxter’s idea of a Q-matrix 
gives an operator construction of the TQFT in terms of q-difference operators, 
so-called quantum D-modules. That this operator construction of a TQFT is 
identical to the one obtained via the Bethe ansatz equations is a non-trivial 
statement: it requires the proof of the existence of an algebra isomorphism 
between the operator version of the TQFT (16.211 . (16.311 and the coordinate ring 
presentation in terms of the Bethe ansatz equation as well as the completeness 
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of the Bethe ansatz; see [H Section 7 and 8]. As a result one obtains an explicit 
algorithm which allows one to compute the fusion coefficients, see [14] and na 
Sec. 7.3.1]. We showed in this article that the fusion coefficients satisfy linear 
systems (1^ and (16.61) with the image of the g-boson algebra generators under 
the quantum Backlund transform. 

In the context of classical integrable systems the connection with TQFT 
rests on the associativity condition of the underlying Frobenius algebra, the 
Witten-Dijkgraaf-Verlinde-Verlinde equations; see Dubrovin’s lecture notes [5]- 
In contrast, associativity is built in from the start in the construction of the 
TQFT in [la as the fusion matrices are realised as endomorphisms over the 
Fock space of the g-boson algebra, that is, the product in the TQFT is the 
ordinary matrix product. The crucial issue here is the commutativity of the 
product which follows from the construction of the commutative subalgebra An 
generated by the quantum integrals of motion of the quantised Ablowitz-Ladik 
chain, {T^} and {Q^}. In other words, for fixed particle number k the 2D 
TQFT is the quantised Poisson subalgebra of classical integrals of motion of the 
Ablowitz-Ladik chain, An,k = article we showed that a 

basis of the quantum integrals of motions, the fusion matrices of the TQFT, is 
obtained from the quantum analogue of Backlund transformations. 

7.3 Open questions 

As pointed our earlier in the text, the main physical significance of our investi¬ 
gation of the Q^-operator for the quantum Ablowitz-Ladik chain is the result 
that its adjoint action describes the discrete time evolution (17:71) of the system. 
To obtain the full time evolution one needs in addition to consider the adjoint 
of the (5+-operator. We wish to exploit this fact in a forthcoming publication 
where we discuss the discrete dynamics of the system and connect it to the 
TQFT. 

An interesting question resulting from our discussion is a possible extension 
of our construction of quantum Backlund transformations to the continuum 
limit, where one obtains the quantum nonlinear Schrodinger model, as well 
as to different boundary conditions; see e.g. [3T], [51] for a discussion of the 
g-boson system with more exotic boundaries. We hope to address these and 
related questions in future work. 
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